Gluons in glueballs: Spin or helicity? 
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In the last decade, lattice QCD has been able to compute the low-lying glueball spectrum with 
accuracy. Like other effective approaches of QCD, potential models still have difficulties to cope 
with gluonic hadrons. Assuming that glueballs are bound states of valence gluons with zero current 
mass, it is readily understood that the use of a potential model, intrinsically non covariant, could 
be problematic in this case. The main challenge for this kind of model is actually to find a way 
to introduce properly the more relevant degree of freedom of the gluon: spin or helicity. In this 
work, we use the so-called helicity formalism of Jacob and Wick to describe two-gluon glueballs. 
We show in particular that this helicity formalism exactly reproduces the J PC numbers which are 
observed in lattice QCD when the constituent gluons have a helicity- 1, without introducing extra 
states as it is the case in most of the potential models. These extra states appear when gluons are 
seen as spin-1 particles. Using a simple spinless Salpeter model with Cornell potential within the 
00 ■ helicity formalism, we obtain a glueball mass spectrum which is in good agreement with lattice QCD 

predictions for helicity-1 gluons provided instanton-induced interactions are taken into account. 
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As quantum chromodynamics (QCD) is built on the nonabelian SU(3)-color group, it allows for purely gluonic 
bound states called glueballs. The structure and properties of these pure glue states is nowadays far from being 
i' completely understood and deserves much interest on both experimental and theoretical sides. 
JL , On the one hand, some experimental glueball candidates are currently known. Most of them are scalar, such as the 
• ao(980), /o(980), /o(1500), /o(1710), . . .but no definitive conclusions can be drawn concerning the nature of these 
states. Indeed, it is often pointed out that the lightest glueballs are probably strongly mixed with other hadrons like 
l— ~~ '■ mesons and tetraquarks for example. Many details concerning the identification of experimental states can be found 
| in the recent report 

On the other hand, pure gauge QCD has been investigated by lattice QCD for many years, leading to a well 
established glueball spectrum below 4 GeV 0,3 ■ Numerous effective QCD models have also been applied to study 
the glueballs. One can quote QCD sum rules 0,01, AdS/QCD correspondence H, Coulomb gauge QCD @, Q, and 
potential models. Pioneering works in this last field were presented in Refs. (9l. ITo|. In both works, glueballs are seen 
as bound states of at least two valence gluons, but the properties of these gluons (mass and spin) are very different. 
In the first reference, valence gluons are assumed to be helicity- 1 particles. It means that their spin has only two 
projections (±1). In the second one however, they are seen as massive particles with spin-1 (the zero projection is 
00 , also allowed). 

A remark should here be done concerning the number of constituent gluons in a given glueball. It appears in lattice 
J> , QCD that the lowest-lying glueballs are the C = + ones. As a bound state of two gluons can only have C = +, it 
is rather natural to assume that the lightest glueballs are mainly two-gluon states (the more constituent gluons are 
present, the more the glueball should be heavy). This picture, that we adopt in the present work, is widely accepted 
in models with constituent gluons. Moreover, it is interesting to mention some results of the Coulomb gauge study 
of Ref . Q . In this approach, a Fock space expansion of glueball states in terms of quasigluons can be performed, and 
it appears that the influence of the three- and four-gluon components on the low-lying C = + glueballs is negligible: 
The two-gluon component is dominant as intuitively expected. But actually, the relevance of using a potential model 
to describe a glueball is still controversial. Let us begin by the problem of the gluon mass. As we already mentioned, 
there are works, in the spirit of Ref. 9], arguing that a valence gluon is a massless particle, which gains a constituent 
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mass /i, either constant 0, or state-dependent = \ yffPj [HI, [13, 13 ■ Relativistic spin-dependent corrections are 

then developed in powers of l//i 2 . In this picture, the valence gluon is a posteriori massive, because it is confined into 
a glueball. Let us note that, more generally, both quarks and gluons can gain a constituent mass from renormalization 
theory (in the Coulomb gauge approach of Ref. [l4j, massless gluons gain a constituent mass of about 0.7 GeV at 
zero momentum). But, other studies keep the assumptions of Ref. 11011 and state that a valence gluon has to be a 
priori considered as massive, that is with a non zero current mass jiol. [l5| . The underlying idea is roughly that 
the nonperturbative effects of QCD cause a mass term to appea r in the gluon propagator. This fixed gluon mass is 
typically assumed to be around m g = 0.5 ± 0.2 GeV [H, fl?i Il8j. and the relativistic corrections are then expanded 
in powers of l/m 2 g . The problem of the gluon spin is obviously linked to its mass. If a valence gluon is a priori 
massive, then it is a spin-1 particle. But, if it is massless, what is the correct internal degree of freedom? The most 
obvious answer is that it has a helicity-1. But, the spin corrections in potential models appear at the order l//^ 2 , at 
a level where the gluons have a dynamical mass, and thus perhaps a spin-1. In the present work, we will only focus 
on the case where the valence gluons are massless, that is have a vanishing current mass. We thus need a formalism 
which allows us to deal with both helicity and spin degrees of freedom, and to build quantum states with the correct 
symmetry following the degree of freedom which is chosen. 

In potential approaches, hadrons are generally described by | 2S+1 Lj) quantum states (in spectroscopic notation) 

which are simultaneously eigenstates of J 2 , L 2 and S 2 , such that J = L + S is the total spin. But actually, J is 
the only relevant angular momentum labeling a hadronic state (together with the parity and the charge conjugation). 
In this picture, it is assumed that a J PC state is a linear combination of the allowed | 2S+1 L,/) states leading to the 
desired value of the total spin. How to build such a general state? We propose here to use the helicity formalism, 
developed by Jacob and Wick in Ref. [l9( to describe scattering in two-body systems. The crucial feature of this very 
general formalism is that it remains valid for massless particles like gluons. But, as we will show through this paper, 
it provides a powerful tool to build J PC states in terms of the usual | 2S+1 Lj) states, and allows to gain considerable 
insight on the glueball models in potential approaches. 

Our paper is organized as follows. Sec. [TT] is a presentation of the helicity formalism. It sums up the key points 
of Ref. [191 ]. Then, in Sec. IIII1 this formalism, which allows for the gluons to have either a helicity or a spin degree 
of freedom, is applied to two-gluon glueballs. The case of massless gluons has already been studied in Ref. [§]. We 
present here a more detailed study with a proper treatment of the relativistic kinematics, and we formulate the glueball 
helicity states in a way that is more convenient to further apply to potential models. Moreover, the glueball spectrum 
is now far better known than at the time of Ref. @ thanks to lattice QCD calculations @, H[. That is why it is of 
interest to reconsider the description of glueballs with helicity states. To this aim, we introduce in Sec. IIVI a simple 
potential model based on a spinless Salpeter Hamiltonian with a Cornell potential. Instanton-induced forces are also 
included. Then, we show in Sec. fVl that the model we introduced, supplemented by the helicity state formalism for 
helicity-1 gluons, leads to a rather good agreement with the lattice QCD spectrum. Finally, we draw some conclusions 
in Sec. ED 



II. HELICITY FORMALISM 

A. General considerations 

Let |V'fe,A) be the state of a particle of momentum k = kl z and helicity A. For a particle of spin s and mass m, one 
can have 

A = -s,-s + l,...,+s (m t^O), (1) 
±s (to = 0). (2) 

If we define a general rotation as R(a,P 7 "f) = exp(— iaJ z ) exp(— i(3J v ) exp(— i'jJz), with a,j3, 7 the Euler angles and 
J the angular momentum operators for the considered particle, then 

|^ P ,A> = R(<t>, 9, -tj>) |^,a) = e l ^i?(0, 6, 0) |V fe ,A) (3) 

is the state of a particle whose helicity is A and whose momentum p (|p| = k) has the arbitrary polar angles (9, cj>). 

In the reference frame where p\ = —pi = p, a two-particle state \ip P ,\ 1 ,x 2 ) can be built from and \i/> p> \ 2 ), 

both given by Eq. It reads 

\iPp,Mm) = \ 1 Pp,>>i)®\Xpm)> ( 4 ) 



with 

lx P ,A 2 ) = (-l) S2 - A2 e-^ 2) |^p,A 2 >- (5) 

The rotation along the y axis ensures that p\ = —p 2 , while the phase factor (— l) s 2->2 j s suc \y that |xo.a) = |"0o,-a) 
as intuitively expected. 

As it is defined in Eq. (0]) , the two-particle state is not an eigenstate of the square total angular momentum J 2 and 
of its projection M(= J z ). However, the state 



\J,M;Xi,X a 



2J+1 



4tt 



1/2 r 2ir 

d<f> I d0smev(i Xi _ X2 (^e : -cf>)R(^e : -cf>)\^ pMM ) (6) 



is, by construction, an eigenstate of J 2 and J z (more details can be found in Ref. [2 (J . Chapter 7] for example). Indeed, 
the Wigner D-functions 

& 7) = ^ lMa d J M, M > (/?) e" 4M ' 7 (7) 

enforce a particular value for the total angular momentum. Their explicit forms can be found for example in Ref. (2ll . 
Chapter 4]. The states given by Eq. ([HJ are orthonormalized by definition, i.e. (J' , M'; X[, X' 2 \ J,M;X%,X 2 ) = 
Sj.j'Sm.m'S^.x^S^,^, and they describe a general two-particle system in the rest frame. The use of helicity de- 
grees of freedom rather than the spin ones allows to deal with massless particles too. This feature will obviously be 
useful in the description of glueballs. Let us note that we did not write explicitly the dependence in p of the helicity 
states in order to simplify the notation. 

B. Symmetries 

The states © are defined so that they satisfy 

J 2 \J,M;Xi,X 2 ) = J(J+1)\J,M;X U X 2 ), (8) 
J z \J,M;Xi,X 2 ) = M\J,M;Xi,X 2 ). (9) 

Moreover, the usual rules concerning the coupling of two angular momenta lead to the constraint 

J>|Ai-A 2 |. (10) 

It can also be shown that the | J, M; X\, X 2 ) states have the following behavior under the parity P: 

P | J, M; Ai, A 2 ) = mmhiy-* 1 -* 2 I J, M; -A x , -A 2 ) , (11) 

where rji is the intrinsic parity of particle i. A physical state is asked not only to be an eigenstate of the total angular 
momentum operators but also of the parity. Such a requirement is fulfilled by the following linear combinations 

\H±; J P ;A 1 ,A 2 ) = ■±={\J,M;\ 1 ,\ 2 )±\J,M;-\ 1 ,-\ 2 )}, A x or A 2 ^ 0, (12) 

\N;J P ) = | J, M; 0,0), (13) 

for which P\H±; J p ;X 1 ,X 2 ) = P \H±; J p ; Ai, A 2 >, with P = ±T ll i l2 (~l) J - Sl - S2 . In the latter, the \H±; J p ; X x , X 2 ) 
and | TV; J p ) states will be referred as helicity states. 

When the two particles are identical (mi = m 2 = m, s± = s 2 = s), it is relevant to study the action of the 
permutation operator P\ 2 . One finds 

[1 + (-l) 2s P 12 ] \ J,M;Xi,X 2 ) = |J,M;A 1 ,A 2 ) + (-l) J |J,M;A 2 ,A 1 ), (14) 

where the operator [l + (— l) 2s Pi 2 ] = S is nothing else than a projector on the symmetric (s integer) or antisymmetric 
(s half-integer) part of the helicity state. Consequently, in the special case of identical particles, the helicity states 
should also be eigenstates of S. By inspection of relations (fTT)f and (fT4"|) . it can be seen that the states 

\H±, P ;J P ) = -^={\H ± -J p -X 1 ,X 2 )+p\H ± ;J p ;X 2 ,X 1 )}, Aj or A 2 £ 0, (15) 
\N-J P ) = | J, M; 0,0), (16) 



with p = ±1, are eigenstates of S. More precisely, 

S\H ± . p -J p ) = [l + p(-l) J ] \H ±tP ; J p ), (17a) 
S\N-J P ) = [l + (-l) J ] \N; J p ). (17b) 

In the symmetric case, we thus observe the emergence of selection rules following the value of J and p. 

C. Wave functions 

It is of great phenomenological interest to be able to express a given helicity state in terms of states of given 
orbital angular momentum L and intrinsic spin S. Indeed, although the total spin J is the only relevant angular 
momentum of the system, especially when one deals with relativistic bound states, most of the Hamiltonian-based 
effective approaches of QCD involve central potentials with relativistic spin corrections. Such Hamiltonians thus act 
on nonrelativistic | 2S+1 Lj) states (in spectroscopic notation) rather than on helicity states. It is actually proved in 
Ref. [lj| that the following decomposition holds: 

1/2 

(L ) S;0,A 1 -A 2 |J,Ai-A2}<* 1 ,a2;Ai,-A2|£,Ai-A 2 ) \ 2S+1 L,,) , (18) 

where we also impose the normalization condition 

( 2S ' +1 L> J ,\ 2S+1 L J )=6 L , L ,5s,s'6 JJ ,. (19) 

The sum (fTg| involves all the {L, S} couples such that S = si+s 2 and L + S = J. The symbols (a, b; c, d\e,f) denote 
the well-known Clebsch-Gordan coefficients. 

More explicitly, one can write the | 2S+1 ij) states as 

(20) 

where the radial, angular, and spin wave functions are explicitly written. It is worth noting that the shape of the radial 
wave function can only depend on the total angular momentum J in order not to destroy the symmetry properties of 
the helicity states. At this stage, one could wonder if the helicity formalism, which has been proposed as a powerful 
way of studying relativistic scattering problems, can be applied to describe bound states. Actually, the construction 
of the helicity states that was presented here is purely geometrical. The angular parts of the states are built in order 
to have the desired properties, but the radial part of each state is arbitrary (excepted that it can only depend on 
J). Consequently, only the dynamics of the system will fix this radial part: Spherical waves will be obtained for 
scattering states, and bound state radial wave functions otherwise. But, the construction we presented up to now, 
being Hamiltonian-indcpendent, is valid in both cases. 

III. HELICITY STATES FOR TWO-GLUON GLUEBALLS 
A. Gluons with helicity 

We now particularize the formalism to the special case of a system made of two gluons with helicity- 1. Then, 
mi = f«2 = 0, Si = «2 = 1, Xi = ±1, r\\r\2 = 1- Actually, the results that we will obtain in this section are 
formally identical to the case of a state made of two photons. The total color wave function is indeed assumed to be 
a singlet one, which is totally symmetric, and does not explicitly appear in the computations. Taking into account 
the symmetrization of the helicity sates, one finds that there are four allowed helicity states, namely 

\S±;J p ) = \H ±tl ;J P ) Xa=Xl , \D + ;J p ) = \H +X ,J p ) X2= _ Xi , \D.;J p ) = \H.,. i; J p ) X2= _ Xi . (21) 

But, the selection rules (fTT]) together with Eqs. (flTJ|) and (TTT|) impose particular values for the total angular mo- 
mentum and parity of these four states. It can be checked that one can only obtain the following states 



|J,M;A 1 ,A 2 )=^ 



r c 



2L + 1 



2J+1 



\ 2S+1 Lj) 



Rj(r))® \Y L (f))<g>\s 1 ,s 2 ) i 



(22) 



The S- and D-labels stand for helicity-singlet and -doublet respectively. 

It should be noticed that a state made of two gluons in a color singlet has always a positive charge conjugation 
(C = +1). More explicitly, the states (f22|) thus give rise to the following glueball states 



|S+;(2*)+) 
\S-;{2k)~) 

\D + ;(2k + 2) + ) 



0++,2++,^ 

2++, 4++,.. 
3++, 5++,.. 



(23a) 
(23b) 
(23c) 
(23d) 



It is readily observed that only the \S±; (2fc) + ) states can lead to J = 0, while the \D±) states always have J > 2. 
Obviously, no J = 1 state is present: Only the \D-) states can generate an odd- J, but J is at least 3 in this case. The 
fact that a state made of two photons (or gluons) can never have the value J = 1 is known as Yang's theorem [2^ |. 
but has also been found independently by Landau (23|. Lattice QCD confirm the absence of 1 h and 1 ++ states, at 
least below 4 GeV. It is worth mentioning that glueball states with an even- J and a positive parity can be built either 
from the helicity-singlet or from the helicity-doublet. The important fact is that the helicity states exactly reproduce 
the J PC glueballs which are observed in lattice QCD without the extra states which are usually present in potential 
models. 

The application of the decomposition formula to the helicity states (|2"2"|) gives 



\S+](2k)+) 

\S-;(2k)-) 
\D+;(2k + 2) + ) 

\D_;(2k + 3) + ) 



1/2 



\ l 2k 2k ) - 
k(2k- 1) 



2fc(2fc + 1) 



3(4fc- l)(4fc + 3) 



-,1/2 



(4fe+l)(4fe-l) 



Hk - 2 2k ) + 



2k 



1/2 



3 2fc-l 2fc )- 



4fc + 1_ 

(fc + 2)(2fc + 3) 1 1/2 
(4fc + 3)(4/c + 5) 



2k + 1 



4fc + 1 



''2fc 2 fc+2) + 



1/2 

| 5 2fc 2fe ) 

(fc + l)(2fc + 1) ' 
(4fc + 3)(4fc + l). 

| 3 2fc + l 2fe ), 

6(fc + 2)(2fc + l)l 1/2 
(4fc + 3)(4fc + 7) 



1/2 



5 2fc + 2 2fe ), 



-,1/2 



(24a) 
(24b) 



5 2/c + 2 2fe+2 ) 



(* + l)(2* + l) 



-,1/2 



(4fc + 5)(4/s + 7) 



2k + 5 
4k + 7 



1/2 



5 2fc + 2 2fe+3 ) 



5 2k + 4 2k+2 ), 

2{k+ 1)1 1/2 



4k + 7 



i 2fc + 4 2fc+3 ) 



(24c) 
(24d) 



TABLE I: Matrix elements of a given operator O for the glueball states composed of two helicity-1 gluons. All non diagonal 
elements are vanishing. The minimal J PC values are indicated in parenthesis. 
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S+;J + ) (0 ++ ) 
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J(J + l) + 2 
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-2 


J(J+1) 
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5_; J") (0-+) 
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J(J + l) + 2 
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-2 


J(J+1) 


2 


D+-J+) (2++) 


1 


J{J + 1) - 2 
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-2 


J(J+1) 


-6 


D„-J+) (3++) 


1 


J(J + 1) - 2 


6 


-2 


J(J+1) 


-6 



Thanks to the decompositions J24|) , the matrix elements of various operators are readily computed. They are given 
in Table HI Note that, through the kinetic energy, the L 2 operator controls the glueball mass in a simple Hamiltonian 
with only a central potential. It appears that the matrix elements only depend on the singlet or doublet nature of an 
helicity state. Moreover, the matrix elements between the various even-J ++ states, \S + ; (J > 2) + ) and \D + - (J > 2) + ), 
vanish. The helicity-singlet and helicity-doublet are thus completely decoupled. It was claimed in Ref. [9( that the 
glueball spectrum should be characterized by a tower of degenerate even- J glueball states with positive and negative 
parity (the helicity-singlets). However, recent lattice QCD computations have unambiguously shown that the ++ 
and ^ glueballs are not degenerate, as well as the 2 ++ and 2 ^ ones @, Q. The difference between the 2 ++ and 
2 ^ states can be easily explained because the lightest 2 ++ states should be an helicity doublet, and not an helicity 
singlet as the 2 '". In addition there should be a 2 ++ helicity singlet with the same mass than the 2 ^ glueball, as 
pointed out in Ref. [§]. This has not been detected in lattice QCD. The problem of the degeneracy of the ^ and ++ 



states however, requires a particular Hamiltonian to be elucidated. One could think for example to instanton-induced 
interactions, which are repulsive in the pseudo-scalar channel, and attractive in the scalar channel [B, H3, HE]. Such 
interactions, that we will further comment in Sec. lIVBj will lead to a correct ordering of the glueball states. 



B. Gluons with spin 



It is also possible that the gluons should be considered as spin-1 particles rather than helicity-1 ones. The most 
obvious way to have this situation is to deal with massive gluons. It should be stressed that, in every glueball potential 
model, the relativistic corrections (containing the spin-dependent terms) involve a constituent gluon mass. This is true 
even if the gluon is massless at the dominant order. Consequently, it could be possible, in the framework of effective 
models, to deal with massless gluons at the dominant order but to give them a spin degree of freedom because they 
are massive at the order of the spin-dependent terms. If the gluons have spin-1 rather than helicity-1 the value = 
can be reached. However, one is always dealing with identical bosons such that m\ — 7712 and 771 = 772. Consequently, 
as it can be deduced from relations (jTTJ) , the helicity states describing a glueball made of two gluons with spin are 
the four states (|2"4"|) supplemented by 

(25a) 
(25b) 
(25c) 
(25d) 
(25e) 

Five additional states appear because of the allowed zero value for the helicity. First, we can point out the apparition 
of a family of glueballs with odd J and negative parity. Such states are actually not observed in lattice QCD. The 
decomposition formula (fT5|) leads to \H+ t —\\ (2k + _ A =1 = \ 3 Jj)> i- e - a pure | 2S+1 ij) state. 
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It is also worth mentioning that a 1 + 
one finds that 



glueball, forbidden with helicity-1 gluons, is now allowed. Applying Eq. (I18D . 



L ;(2fc + l)+) Ai= 



0,A 2 = 1 



2k 



4k + 3 



1/2 



'2k 



2k+l 



2k + 3 
4k + 3 



1/2 



5 2fc + 2 2fc+ i), 



(26) 



that is a decomposition which is very similar to the helicity-doublet j24d|, which generates the same quantum 
numbers for J > 3. Both states can actually couple to each other through the orbital angular momentum: 
(H— —i; J + \ L 2 |.D_; J + ) = 2y/ (J — 1)(J + 2). Let us further investigate this point by considering the dynamics 
of the system through its Hamiltonian H. As the |-D_; J + ) and \H— —\\ J + ) states are coupled, the physical states 
are eigenvectors of the Hamiltonian 



n 



(£>_; J+\H |D_;J+) 
J+\H\D_;J^ 



(D-;J+\H\H-,-i;J+) 
H-,-i;J+\H\H--i;J+ 

/2s 1 



(27) 



where these four matrix elements can be expressed as linear combinations of ^ 2S +1 L'j H | 2S+1 _L,/) thanks to rela- 
tions P4"|) and |^|) . Then, it can be shown that the eigenstates of Hamiltonian (f2"Tf are the pure | 5 J— lj) and 
| 5 J + 1,/^ states. This result is obtained under the assumption that no tensor force is present. This is always the case 
at the dominant order, in particular when one deals with central potentials with relativistic corrections. 



In the same way, assuming that ( 2S+L L' J H\^ ti+1 Lj) cx 5s'.s5l',l, one can check that the coupling between 
the \S- 4 , J + ) and J + ) Xi= q a 2 =i s t & t es f° r J > 2 leads to the conclusion that | 3 J — lj) and | 3 J + 1,/) are the 

quantum states that have to be considered. Finally, the \N; J + ) and |-H+,i; J + ) Ai=0 x ± and helicity states are 
coupled to each other and to the \S+; J + ) and \D+; J + ) states, so that we have checked that the physical quantum 
states are I 1 Jj), | 5 J-2j), | 5 Jj), and | 5 J + 2j). 

In conclusion, when the dynamics of the system is included, the nine possible helicity states reduce to the nine 
possible \ 2S+1 Lj) states that are usually used in potential models of glueballs. Actually, it is rather logical that the 
helicity formalism reduce to a usual LS'-basis when particles with spin are considered since all the spin projections 
are allowed. Let us remark that for J = 0, the two physical states are given by ^So) and | 2 -Do)- As the v state is 
a pure \ 3 Po) state, the nondegeneracy of the scalar and pseudoscalar glueballs is de facto explained when the valence 
gluons have a spin degree of freedom. 



IV. POTENTIAL MODEL 



A. 



Main Hamiltonian 



The construction of two-gluon helicity states presented in Sec. Mil is based on purely kinematical arguments. In 
order to compute a glueball mass spectrum, it is necessary to use a particular Hamiltonian which will contain the 
dynamics of the system. The simplest way of modeling a two-gluon glueball is to use a two-body spinless Salpeter 
Hamiltonian with a Cornell potential, that is 



The kinetic part is the kinetic energy of two spinless massless particles, i.e. the valence gluons for which the spin is 
neglected at the dominant order. But, the spin symmetry will be taken into account by the use of helicity states, even 
if Hamiltonian (|28p is spin-independent. 

The potential term has a Cornell shape, that is a linear-plus-Coulomb form. The linear confining term can be 
seen as the static energy of a flux tube linking the two gluons. The string tension a g can be related to the string 
tension of a mesonic flux tube, denoted as a, by a scaling law a g = Co. A typical value for a is about 0.2 GeV 2 , 
and two values of C can be found in the literature: either C = 9/4 (Casimir scaling), or 3/2 (square root of Casimir 
scaling). While the 3/2 factor is commonly found in bag model-inspired approaches [26], the Casimir scaling seems 
to be favored by more recent effective approaches and by lattice computations [13] . We will here assume the Casimir 
scaling hypothesis, that is C = 9/4. Beside the long range linear potential, the Coulomb term comes from short range 
interactions: It is the lowest order approximation of the one gluon exchange diagram between two gluons. a s is then 
an effective strong coupling constant smaller than 1 and the factor 3 is the color factor associated with a gluon pair 
in a color singlet. It was shown in Ref. [28| that, starting from the ++ glueball mass and wave function as computed 
in lattice QCD, the inverse problem can be solved, and the equivalent Hamiltonian is compatible with the form (|28p . 
This validates such an Hamiltonian description, at least in the case where the valence gluons have a spin, since in 
Ref. [28[ it is assumed that the scalar glueball is a ^Sq) state. 

Relativistic corrections to the Cornell potential can also be computed. For example, in the flux tube model, 
corrections to the linear potential appear as a dynamical term proportional to a g L 2 and a spin-orbit term proportional 
to a g L ■ S [1^] . A nonperturbative retardation term has also been proposed [3(| ■ Moreover, relativistic corrections 
to the Coulomb term can be computed from the QCD Feynman diagrams involving two gluons at tree level. Their 
complete expression can be found in Ref. [3l| . and involves the usual contact (spin-spin), spin-orbit, and tensor 
interactions. For our purpose, it is sufficient just to list the global structure of all these additional terms. All the 
matrix elements appearing in these tree level relativistic corrections can be found in Tables [J for helicity states. 
Finally, it is worth mentioning that the first relativistic corrections are, in this formalism, proportional to l//i 2 , /i 

being the dynamical mass jj, = ( / gained by the valence gluons because of confinement. 



It is known that, in the light meson sector, nonperturbative contributions due to instantons comes into play. Roughly 
speaking, instantons are classical solutions of the Euclidean equations of motion of QCD, which provide informations 
on the nontrivial vacuum structure of this theory. We refer the reader to Ref. [32j for a review on instantons in 
QCD. It has been shown that instantons-induced forces exist between the quark and the antiquark in a light meson. 
Such forces can be included in potential quark models as an isospin-dependent contact term which is nonzero in the 
pseudoscalar channel (0 *") only 33]. In particular, the strong attractive nature of the instanton-induced contribution 
in the case of the pion is able to explain the particularly low value of its mass, without lowering the masses of the 
other mesons. 

If the instanton-induced forces are rather well understood in the meson sector of QCD, the situation is not so clear 
for glueballs. It has firstly been shown in Ref. |2J] that instantons induce a strong attractive force in the scalar glueball 
channel and a repulsive force in the pseudoscalar channel. In the tensor (2 ++ ) channel, these forces vanish; moreover, 
instantons are not expected to play any role in the other channels. From the results of a more recent study [5] , it is 
tempting to assume that the instanton-induced forces in the scalar and pseudoscalar channels are of equal magnitude 
but of opposite sign. This could be a consequence of the self-duality of the instanton's field strength and a general 
characteristic of the instanton contribution in all hadrons [25[ . From this discussion, and although its exact form has 
not been computed yet, we can propose the following ansatz for the instanton-induced contribution: 




(28) 




B. Instanton-induced forces 



AH! = -PI 6 



(29) 



Such a term only contributes for J — and depends on the parity P. Its magnitude is related to the unknown 
parameter X, that we assume to be positive and constant in first approximation. 

Why is it so interesting to study the influence of instantons on our model? The problem actually comes from the 
important mass splitting between the ++ and h glueballs that is observed in lattice QCD. If the valence gluons are 
spin-1 particles, instanton-induced interactions are not needed to explain this nondegeneracy, as it has been argued in 
Sec. lIIIBl But, if the valence gluons are helicity-1 particles, then an additional mechanism is required to lift the mass 
degeneracy of the scalar and pseudoscalar glueballs. As it can be seen in Table HI no correction involving the usually 
used operators will be able to do that, since their matrix elements are identical for the ++ and h glueballs. That 
is why instanton-induced forces are particularly interesting. First, they act in the correct way, increasing the 
mass and decreasing the ++ one. Second, they have already proved to be very useful in the meson case, and, since 
we know from Refs. 0, [HJ that instantons play a role in glueballs too, it would seem more coherent if their effects 
were included in a glueball model also. 



V. MASS SPECTRUM 



A. Parameters 



Before performing explicit computations, it is necessary to fix the different parameters of our model. As we 
already said, we focus on valence gluons with vanishing current mass and we assume the Casimir scaling hypothesis: 
(9/4)er. We set a — 0.185 GeV 2 for the mesonic string tension. This value is located in a rather standard 



a 



<j 



interval: a S [0.17,0.2] GeV 2 is commonly found in the literature. Moreover, this particular value has already given 
very good results in a previous computation of quarkonium mass spectra in the flux tube model [l3l ] . Two models will 
be proposed following that the valence gluons are assumed to be spin-1 (Model A) or helicity -1 ( Model B) particles. 



In Model A, the ++ ground state is a L = S = one. In this case, we have shown in Ref. 28] that a s — 0.2 was 
compatible with the current lattice QCD data. No instanton-induced interaction is needed in this case since the scalar 
and pseudoscalar glueballs are de facto nondegenerate. We thus set X = 0. In Model B, the situation is more similar 
to the Coulomb gauge approach of Ref. [a], since gluons with helicity are used. We will take a s = 0.45, a value close 
to the one of Refs. 0, In this case, instanton-induced interactions are required, and we set X = 0.45 GeV in order 
to reproduce at best the ++ and ^ masses which are computed in lattice QCD [1,0]. The parameters appearing 
in both models are summed up in Table ITT1 



TABLE II: Parameters used in our computations. 





Model A 


Model B 










a (GeV 2 ) 


0.185 


0.185 


a s 


0.200 


0.450 


1 (GeV) 





0.450 



The glueball mass spectrum now remains to be numerically computed from the central spin-independent Hamilto- 
nian Hq. Only the radial wave function will be affected by the Hamiltonian of the system, the spin and angular parts 
being fixed thanks to the helicity formalism. We will use the Lagrange mesh method to compute the matrix elements 
of Hq. This method allows a very simple and accurate treatment of semirelativistic Hamiltonian of the form (|28p . 
We refer the reader to Refs. [34], [35j for more informations about the Lagrange mesh method. Knowing the matrix 
representation of Hq, the eigenequation 

H Q \j PC )=M \j PC ) (30) 

can be solved, and the total mass is given by 

M = M -PXSj, . (31) 



B. Results 



A glueball mass spectrum can now be computed within the framework of our spinless Salpeter model with the 
helicity formalism. The results will be compared to some lattice QCD predictions concerning the glueball masses (see 



TABLE III: Available data for C = + glueball masses from various lattice QCD models and Coulomb gauge QCD (CGQCD), 
compared with the results of Models A and B. The glueball mass is given and the corresponding spin/helicity state is detailed 
in both cases. Parameters of Table [TT] are used, and all masses are given in GeV. 
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Table Hill and Fig. [T]). The parameters of our models are fitted on data taken from Ref. [3| completed by data from 
Refs. [2l l36| . Results from Ref. are also given. The predictions of this reference are quite different from the results 
coming from the compilation of Refs. 0, Q . The ground states of Ref. [37| have lower masses, and more excited states 
have been computed. 

It is also interesting to compare our results with the Coulomb gauge approach of Ref. [8J. In this last reference, 
helicity-1 gluons are considered and encoded in the Coulomb gauge formalism. The main features of the mass spectrum 
of Ref. @ should then be similar to our Model B. This is roughly the case but with a serious exception: The mass 
gap between ++ and ^ states is about 200 MeV, which is far less than the value predicted by lattice calculations. 
But no instanton-induced interaction is explicitly taken into account in this work. 

A detailed glueball spectrum is given in Table ITTT1 for Models A and B. We computed the masses of more states 
than those which are currently observed in lattice QCD. Some of them have a mass greater than 4 GeV; glueball 
spectrum in lattice QCD is poorly known above this energy range. However, there are higher ++ and 2 ++ states 
that lie under this limit with both sets of parameters. Some of them are seen in Ref. [37} but not in Refs. @, It 
should be interesting to know whether future lattice computations will confirm or not the existence of these states. 
We also point out again that no J = 1 state is present at low energy as expected from lattice QCD with helicity-1 
gluons. 

Let us begin by a discussion of the results obtained with Model A. In this case, the scalar and pseudoscalar glucballs 
are compatible with lattice QCD without invoking instanton-induced interactions. But, the situation gets clearly worse 
for higher J. First, J — I states are present, which are not observed in lattice QCD. In particular, the rather light 
1 h glueball seems to be a serious flaw of Model A since the — + channel is rather well known from lattice QCD in 
this energy range, and no such state has been seen. Then, the lightest 2 ++ state is degenerate with the ++ state, 
and it should not be the case. Actually, nearly every state with Model A does not lie within the error bars of lattice 
QCD, suggesting that this model should be discarded. 

We turn now our attention to Model B. In this case, an instanton-induced term is needed, otherwise the ++ and 
^ glueball would have the same mass. The value X — 0.450 GeV is of the same order than the typical magnitude 
of instanton-induced effects in mesons [32| ]. Globally, the rest of the spectrum is in agreement with lattice QCD, and 




FIG. 1: Comparison between the lattice QCD data concerning C — + glueballs (crosses), the Coulomb gauge results (trian- 
gles) [U], and our Model B (circles). Masses are given in GeV. All lattice data come from Refs. @, Oil except data for 4 ++ and 
6 ++ states (37J (see 2nd and 3rd columns of Table ITTTjl . 



TABLE IV: Comparison of C = + glueball mass ratios, normalized to the lightest ++ state, between the lattice QCD 
data [3, Hli and our Model B (the corresponding helicity state is given). Parameters of Table HT1 are used. The lattice mass 
ratios followed by Q are taken from Table VIII of this last reference. The other ones, that were not given in Ref. Q, have been 
computed with the more recent data of Ref. |3j. 
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the agreement is far better than with Model A. Our (2fc + 2) ++ states, although being roughly compatible with lattice 
QCD, lie in the upper part of the errors bars, while the 3 ++ state is too light. A modification of the parameters a or 
a s would shift the whole spectrum in the same direction. So an improvement of the 3 ++ mass would spoil the rest 
of the spectrum. Spin-dependent interactions at tree level (see TableHJ are the same for \D+; J + ) and \D—; J + )- So, 
these interactions would shift the 3 ++ and 2 ++ masses in the same direction. We suggest that additional mechanism 
such as singlet-doublet mixing - corrections beyond tree level - could cure this mass problem. Let us note that, in 
Ref. 0, even though the 2 ++ is located inside the lattice error bars, the 3 ++ is also below. Model B is thus rather 
convincing, particularly because the number of states is drastically decreased with helicity- 1 gluons: The few states 
that are observed in lattice QCD are the only one that can be built, without extra states as it is the case in Model A. 
The mass spectrum of Model B has been plotted in Fig. [T] and is compared to lattice QCD and Coulomb gauge data. 

The errors of lattice data on absolute glueball masses are quite large. This is due the the uncertainty in setting the 
mass scale. This problem can be corrected by computing mass ratios. We choose to report all masses to the lowest 
++ state as in Ref. @, Table VIII]. This state does not present the lowest statistical uncertainty but it is the ground 



state. The mass ratios for our model B, computed with data from Table IPIII are given in Table HVl and Fig. and 
are compared with the corresponding lattice data [H, Q • The study of mass ratios brings the same conclusions as the 
discussion above about absolute masses. The mass ratios obtained with Model A are not mentioned since this model 
does not bring relevant results. 
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FIG. 2: Comparison of C = + glueball mass ratios, normalized to the lightest ++ state, between the lattice QCD data 
(crosses) @, H|, and our Models B (circles). More details on these data can be found in Tables iLUl and ITVl 



VI. CONCLUSIONS 

We have presented in this work an application of the helicity formalism developed in Ref. pi! for bound states of two 
valence gluons. We have shown that this helicity formalism only leads to the quantum numbers which are observed 
in lattice QCD. This feature is a considerable improvement of usual potential models in which lots of extra glueball 
states are obtained. In particular, the well-known absence of J = 1 states has been shown to be a consequence of the 
symmctrization in the helicity formalism for helicity-1 gluons. J = 1 states are however allowed for spin-1 gluons. 

We have then developed a simple potential model, relying on a spinless Salpeter Hamiltonian with a Cornell 
potential and massless gluons. However, we stress that our framework, based on the helicity formalism, can be 
applied to any existing potential model provided that the helicity states are used to compute the different matrix 
elements. An extra physical mechanisms has been considered: instanton-induced interactions, for which we have 
proposed a phenomenological term taking into account their main properties. Using the Lagrange mesh method, 
the glueball spectrum coming from our Hamiltonian model can be computed within the helicity formalism. Two 
possibilities were taken into account: Either the valence gluons are spin-1 or helicity-1 particles. We have computed 
the masses of different states by using standard values for the string tension and the strong coupling constant. The 
instanton parameter has been fitted so that it leads to an optimal agreement with lattice QCD data in the scalar 
and pseudoscalar channels. It appears that, if gluons are spin-1 particles, no instanton-induced term is needed in the 
scalar and pseudoscalar channels, but the rest of the spectrum is not in agreement with lattice QCD. If helicity-1 
gluons are assumed, an instanton term has to be added since the ±+ glueballs would be nondegenerate otherwise. 
Its value is similar to the one encountered in the meson sector. In this case, the glueball mass spectrum is in good 
agreement with lattice QCD, both qualitatively (no extra J p+ states are obtained) and quantitatively. 

By comparing the results of both pictures with lattice QCD, it appears that modeling the currently known C = + 
glueballs by a bound state of two massless valence gluons with helicity-1 seems to be far more relevant, justifying 
a posteriori the assumption that two-gluon states dominate in this sector. In particular, the necessity of adding 
instanton-induced forces should not be seen as a flaw of the model, but rather as a way to be more coherent with 
other studies showing that instantons contribute in glueballs. 



In conclusion, the helicity formalism appears to be a very promising way to improve potential models of hadrons 
containing valence gluons, because it takes into account correctly the relativistic character of these particles. It is 
relevant to assume that such valence gluons are massless particles with helicity-1, as we argued in this paper. For what 
concerns the glueball mass spectrum, it is remarkable how the addition of the helicity formalism to a simple potential 
model as the one we developed here leads to such a nice agreement with lattice QCD. In future works, we plan to 
present a more accurate glueball model, still within the helicity formalism but including relativistic corrections. 



Acknowledgments 

F. B. and C. S. thank the F.R.S.-FNRS for financial support. V. M. thanks the IISN for financial support. The 
authors are grateful to N. Boulanger and T. Schaefer for valuable discussions and advice about the present work. 



[1] E. Klempt and A. Zaitsev, Phys. Rept. 454, 1 (2007) ^^1^:0708.4016, and references therein 



[2] C. J. Morningstar and M. Peardon, Phys. Rev. D 60, 034509 (1999) hepdat/9901004 

[3] Y. Chen et al, Phys. Rev. D 73, 014516 (2006) hep-lat/05 10074 . 

[4] S. Narison, Nucl. Phys. Proc. Suppl. 96, 244 (2001) hep-ph/0012235 

[5] H. Forkel, Phys. Rev. D 71, 054008 (2005) [hep~p h/0312049|T^ 

[6] H. Boschi-Filho, N. R. F. Braga, and H. L. Carrion, Phys. Rev. D 73, 047901 (2006) hep-th/0507063 



[7] A. Szczepaniak, E. S. Swanson, C. R. Ji, and S. R. Cotanch, Phys. Rev. Lett. 76, 2011 (1996) hep-ph/95 11422 . 

[8] A. P. Szczepaniak and E. S. Swanson, Phys. Lett. B 577, 61 (2003) [hep^p h/0308268] . 

[9] T. Barnes, Z. Phys. C 10, 275 (1981). 

[10] J. M. Cornwall and A. Soni, Phys. Lett. B 120, 431 (1983). 

[11] A. B. Kaidalov and Yu. A. Simonov, Phys. Lett. B 636, 101 (2006 ) |hep-ph/051215ll . 

[12] F. Brau and C. Semay, Phys. Rev. D 70, 0140 17 (2004) [hep-ph/0412173| . 

[13] F. Buisseret, Phys. Rev. C 76, 025206 (2007) [arXiv:07 05.0916 . 



[14] F. J. Llanes-Estrada, P. Bicudo, and S. R. Cotanch, Phys. Rev. Lett. 96, 081601 (2006) hep-ph/0507205 . 
[15] W.-S. Hou and A. Soni Phys. Rev. D 29, 101 (1984); W.-S. Hou, C.-S. Luo, and G.-G. Wong, Phys. Rev. D 64, 014028 
(2001) [hep-ph/OlOlTIo] . 

[16] C . W. Bernard, Phys. Lett. B 108, 431 (1982); A. C. Aguilar and J. Papavassiliou, JHEP 0612, 012 (2006) 
[hep-ph/0610040] 



[17] F. D. R. Bonnet, P. O. Bowman, D. B. Leinweber, and A. G. Williams, Phys. Rev. D 62, 051501 (2000) hep-lat/0002020 . 

[18] A. C. Aguilar, A. Mihara, and A. A. Natale, Phys. Rev. D 65, 054011 (2002) [hep-ph/0 109223 . 

[19] M. Jacob and G. C. Wick, Ann. Phys. 7, 404 (1959). 

[20] J. R. Taylor, Scattering theory (John Wiley and Sons, 1972). 

[21] D. A. Varshalovich, A. N. Moskalev, and V. K. Khersonskii, Quantum Theory of Angular Momentum (World Scientific, 

Singapore, 1988). 
[22] C. N. Yang, Phys. Rev. 77, 242 (1950). 

[23] L. D. Landau, Dokl. Akad. Nauk. SSSR 60, 207 (1948); for an english summary of this reference see also H. Fritzsch and 

P. Minkowski, II Nuovo Cimento 30, 401-402 (1975). 
[24] T. Schafer and E. V. Shuryak, Phys. Rev. Lett. 75, 1707 (1995) hep-ph /94103721 . 
[25] N. Kochelev and D. P. Min, Phys. Lett. B 633, 283 (2006) arXiv:hep-ph/0508288 . 

[26] K. Johnson and C. B. Thorn, Phys. Rev. D 13, 1934 (1976); G. Martens, C. Greiner, S. Leupold, and U. Mosel, Phys. 
Rev. D 70, 116010 (2004). 

[27] L. Del Debbio, M. Faber, J. Green site, and S . Olejnik, Phys. Rev. D 53, 5891 (1996) hepdat/9510028 ; G. S. Bali, Phys. 

Rev. D 62, 114503 (2000) [hepTa t/0006022 ; C. Semay, Eur. Phys. J. A 22, 53 (2004) hep-ph/0409105 ; P. Bicudo, M. 

Cardoso, and O. Oliveira. larXiv:0704.2156l 
[28] F. Buisseret and C. Semay, Eur. Phys. J. A 33, 87 (2007) [hep-ph/0611216|. 
[29] F. Buisseret and C. Semay, Phys. Rev. D 76, 017501 (2007) arXiv:0704.1753l. 
[30] F. Buisseret and C. Semay, Phys. Rev. D 72, 114004 (2005) hep-ph /0505168 r . 
[31] V. Mathieu and F. Buisseret, J. Phys. G. 35, 025006 (2008) hep-ph/0702226 . 
[32] T. Schafer and E. V. Shuryak, Rev. Mod. Phys. 70, 323 (1998) hep-ph/9610451 . 
[33] F. Brau and C. Semay, Phys. Rev. D 58, 034015 (1998) [hep-ph/0412179| , and references therein. 
[34] D. Baye and P.-H. Heenen, J. Phys. A 19, 2041 (1986); D. Baye, J. Phys. B 28, 4399 (1995). 

[35] C. Semay, D. Baye, M. He sse, and B. Silve stre-Brac, Phys. Rev. E 64, 016703 (2001); F. Buisseret and C. Semay, Phys. 

Rev. E 71, 026705 (2005) [hep-ph/0409"033] . 
[36] D. Q. Liu and J. M. Wu, Mod. Phys. Lett. A 17, 1419 (2002) 
[37] H. B. Meyer and M. J. Teper Phys. Lett. B 605, 344 (2005) | 

Ph.D. thesis hepdat/0508002 . 



hep-lat/0105019 



iep-ph/0409183 ; H. B. Meyer, Glueball Regge trajectories, 



